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Abstract 

Starting from a representation of the early time evolution of a dynamical system in terms of the polynomial 
expression of some observable c/>(t) as a function of the time variable in some interval < t < T, we investigate 
how to extrapolate/forecast in some optimal stability sense the future evolution of <f>(t) for time t > T. Using 
the functional renormalization of Yukalov and Gluzman, we offer a general classification of the possible regimes 
that can be defined based on the sole knowledge of the coefficients of a second-order polynomial representation 
of the dynamics. In particular, we investigate the conditions for the occurence of finite-time singularities from the 
structure of the time series, and quantify the critical time and the functional nature of the singularity when present. 
We also describe the regimes when a smooth extremum replaces the singularity and determine its position and 
amplitude. This extends previous works by (1) quantifying the stability of the functional renormalization method 
more accurately, (2) introducing new global constraints in terms of moments and (3) going beyond the "mean-field" 
approximation. 
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1 Introduction 



Finite-time singularities in the dynamical equations used to describe natural systems are not always pathologies that 
should be thrown away or ignored but may often betray important useful information on the characteristic properties 
of the real system. Actually, spontaneous singularities in ordinary and partial differential equations are quite common 
and have been found in many well-established models of natural systems, either at special points in space such as in 
the Euler equations of inviscid fluids [Q, in the surface curvature on the free surface of a conducting fluid in an 
electric field in vortex collapse of systems of point vortices in the equations of General Relativity coupled to 
a mass field leading to the formation of black holes [Qj, in models of micro-organisms aggregating to form fruiting 
bodies or in the more prosaic rotating coin (Euler 's disk) [{7|, |||. Some more complex examples are models of 
rupture and material failure [f| [H], 11], earthquakes [12, |3|] and stock market crashes [Jl4|, Op. 



In a recent work [Jig], we have developed theoretical formulas for the prediction of the singular time of systems 
which are a priori known to exhibit a critical behavior, based solely on the knowledge of the early time evolution of 
an observable. From the parameterization of such early time evolution in terms of a low-order polynomial of the time 
variable, the functional renormalization approach introduced by Yukalov and Gluzman allows one to transform this 
polynomial into a function which is asymptotically a power law. The value of the critical time t c , conditioned on the 



assumption that t c exists, can then be determined from the knowledge of the coefficients of the polynomials. Ref.Qlq] 
has tested with success this prediction scheme on a specific example and showed that this approach gives more 
precise and reliable predictions than through the use of the exact real power law model. This is a rather surprising 
and paradoxical observation in contradiction with common wisdom according to which the best and most reliable 



prediction should be obtained from the use of the exact underlying model of the dynamical behavior Q17Q. The reason 
why this is not always true is that our method has shown that approximate solutions can be more stable and more 
reliable when devised so as to maximize the criterion of stability with respect to perturbations in the functional space 



(Ml. 

Here, we extend this work by offering a general classification of the possible regimes that can be defined based 
on the sole knowledge of the coefficients of the polynomial expansion of some observable as a function of the time 
variable. We mostly restrict our analysis to second-order polynomials. As a consequence, the classification can be 
organized in a unique way in terms of the single signed Froude parameter, ratio of the square of the velocity to the 
acceleration. 

Let us assume that the dynamical behavior of a system is sampled to obtain a time series 4>{t) in the interval 
< t < T. The question we address in this work is how to extrapolate/forecast in some optimal sense the future 
evolution of cj)(t) for time t > T. Beyond a simple extrapolation, we ask whether it is possible to detect the germs of a 
finite-time singularity from the structure of the time series, and quantify the critical time and the functional nature of 
the singularity when present. We also aim at classifying the regimes when a smooth extremum replaces the singularity 
and at determining its position and amplitude. 

We shall work in the framework of the functional renormalization method, which constructs the extrapolation for 
t > T from a resummation of the time series represented by a simple model-free polynomial expansion in powers of 
time t, where t is counted since the beginning of the recorded time series. 

2 Summary of the functional renormalization approach to extrapolation of times 
series 

The mathematical foundation of the functional renormalization approach used here can be found in earlier publications 



[]T8[] - [pQ[| , which we summarize briefly. Let perturbation theory (or some fitting procedure) give for the time series <p(t) 



the succession of non-random approximations 4> n (t) where n = 0,1, 2, ... enumerates the order of the increasingly 
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precise approximations. The case of a polynomial expansion will be discussed below, 

n 

<^) = 5> fc i fc , n = 0,1,2,3,... (1) 

A:=0 



This expansion (jlj) is in principle defined for t sufficiently small and based on information about the time series up 
to the time T. The expansion may have no meaning if continued straightforwardly to the region of finite t > T. 
The problem of reconstructing the value of the function in some distant moment of time from the knowledge of 
its asymptotic expansion as t — > is called a renormalization or resummation problem in theoretical physics. An 
analytical tool for the solution of this problem, called algebraic self-similar renormalization, has been developed 



recently J18P-Q3O0 of which we summarize the salient points useful for the present work. 

It is convenient to remove the constant term and consider the time series p(t) = 4>(t) — <po represented by the 
sequence of polynomial approximations Pi(t),i = 1, 2, n, 

Pl (t) = ai t, p 2 (t) =pi(t) + a 2 t 2 , p n (t) = p n -i{t) + a n t n . (2) 

The algebraic self-similar renormalization starts by applying to the approximations (Q) a simple algebraic multiplica- 
tion, thus defining a new sequence, Pi(t, s) = t s pi(t), i = 1,2, n, with s > 0. This transformation increases the 
powers of the series ([I]) and (^). This formal manipulation effectively increase the order of the expansion and provides 
a trick to effectively take into consideration more points from the system trajectory. In the first part of the paper, we 
use the strongest form of this transformation corresponding to formally take the limit s — > oo. This will be shown to 
provide a representation of the renormalized function as an imbedded set of exponentials associated with a singularity 
<ft ~ (*c — t) z > w i m critical exponent z imposed to have an absolute value equal to 1. This limit can thus be considered 
as analogous to a "mean field" regime with the mean field value of the critical exponent z = — 1. This can for instance 
represent the exponent of the derivative of a quantity exhibiting a weak logarithmic divergence. In the second part of 
the paper, we shall relax this limit and determine both t c and the index z self-consistently. 

The second step consists in considering the sequence of transformed approximations Pj(i, s), as a dynamical 
system in the discrete "order time" equal to the order i = 0, 1, n — 1 of the approximation. In order to keep the 
information on the system evolution with the real time t, it is convenient to introduce a new variable tp and define the 
so-called expansion function t(ip, s) from the equation P\(t, s) = ait 1+s = tp, which gives t(<p, s) = (tp/ 'ai) 1 ^ 1+s . 
We then construct the discrete flow in the space of approximations indexed by the "order time" as 

yi(p,s) = Pi(t(ip,s),s) . (3) 

One can then write the equation of evolution in the space of approximations as a function of the discrete "order time" 
in the form of the functional self-similarity relation 

y i+p {p, s) = yi(y p (<p, s), s) . (4) 

Expression (Q) provides the necessary condition for the self-consistency of the cascade of discrete approximations and 
ensures the convergence of Pi's. Expression (Q) is nothing but the application of the Caccioppoli-Banach principle 
for the existence of a stable fixed point (see for instance Chapter XVI of [|3l|]). 

At this stage, the efficiency of the algebraic transformation can be checked by analyzing its stability. This is done 
for the sequence of Pi(t, s)'s by calculating the so-called local multipliers (essentially proportional to the exponential 
of Lyapunov exponents) 

_ dyj(ip.s) 

mi(t,s) = [ — — — \ v=Pl(ttS) . (5) 

When all \rrii(t, s)\ < 1, the convergence of the sequence Pi is guaranteed. To implement concretely the calculations, 
we use the integral form of the self-similarity relation (Q) 

dp 



P»_l Vi(p,s) 



(6) 
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where the cascade velocity is Vi(tp, s) = yi(ip, s) — ]_(<£>, s) and r is the minimal number of steps of the approxi- 
mation procedure needed to reach the fixed point P*(t, s) of the approximation cascade. It is possible to find P*(t, s) 
explicitly and to perform an inverse algebraic transform after which the limit s — > oo is to be taken. This completes 
a first loop of the self-similar renormalization. This procedure can be repeated as many times as it is necessary to 
renormalize all polynomials which appear at the preceding steps. This is the main idea of the self-similar bootstrap 



Q22p- 

Completing this program, we come to the following sequence of self-similar exponential approximants 

Pj(t,Ti t T2,...,T j-x) = aiiexp ^tr x ...exp y-j^-Tj-xtj J j = 2,3..., n . (7) 

and 

0i(*,Tl,T2,-,Ti-l) =Pj(t,T 1 T 2} ...,T j -. 1 ) + a . (8) 

Explicitly, for the three first orders, we obtain 

^(^l) = ai^exp (^ tr i) + a o , (9) 

<#j(*> T i> T 2) = aitexp ^tnexp ^~* T 2^ + a o ; (10) 
<t>l(t,n,T 2 ,T 3 ) = aitexp ^tnexp (j^tr 2 exp ("^s))) + a o • (H) 



Compared to previous works Q18|]-[[30|], we introduce the following innovation to improve on the selection of the 
stable extrapolations (or scenarios). In order to check whether the sequence of 4>j{t, Tx t T 2> ..., t j-x) really converges, 
in addition to calculating the multipliers rrii(t,s) defined by (Q), we analyze the multipliers Mj(t,Tx t T2 t ., ,Tj-x) 
corresponding specifically to <ftj(t, txT2,.-,t j-x)- For this, we construct again an approximation cascade as described 
above and define 



Mj{t,rxT 2r , Tj-x) = [— ^ l=p l(t ,o) • (12) 



Such a definition of multipliers allows us for the first time to compare the convergence of the expansion and of the 
renormalized expansion, making clear the improvements that can be expected a priori from the technique. If the 
values of Mj are systematically smaller than the corresponding rrij in the region of t < T and the overall convergence 
properties are improved, one can expect that the renormalized expressions will work better than the original expansion 
at t > T. 

The final step consists in determining the control parameters tx ; t 2i ... t j-x, by expanding 4>*At, txt 2i ...,t j-x) m 
the vicinity of t = 0, and requiring that this expansion agrees term-by-term with the initial one in terms of the <f>j{t). 
For each j, we obtain j — I self-similar approximants for the sought function (where all control parameters r are now 
known functions of the parameters a, with t\ = 1), 

0*i(M,l-,l) , 

<f>* j2 (t,l,T 2 l,...,l) , 

(13) 

0* j _ 1 (t,l,r 2 ,..,r i _i) , (14) 

which differ according to the number of control functions being used. We can now construct a table of self-similar 
approximants, varying both j and the number of controls. Accordingly, we can define the table of multipliers, varying 
both j and the number of controls. For instance, for j = 4, we have the following table of approximations: 

<j>* 21 {t) = 4>*{t,i) , 
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^ 1 (i) = ^(i,l,l), ^ 2 (i) = ^(i,l,r 2) ), 

0* x (t) = 0*(t, 1,1,1), # 2 (i) = $(t,l,7i,l), (t) = 1,t 2 ,t 3 ) . 

In order to select the scenarios for the extrapolation to the future from the initial time series, one needs to examine 
both the properties of convergence of the sequences of multipliers M, and of approximants <fi*. The best scenario is 
the one which exhibits the best simultaneous convergence of both sequences. The resulting limiting fixed point read 
from the table of approximants should then be taken for the sought extrapolation function. If there are more than one 
limiting points, the sought function should be constructed by taking their average with weights inversely proportional 



to the multipliers (see []30|]). We now illustrate this general methodology on specific examples for which a finite-time 
singularity might exist. 



3 Classification of singular and non-singular behaviors based on second order ex- 
pansions 

3.1 Definitions 

Let us now consider the simplest non-trivial case allowing for the possible existence of a finite-time singularity, 
namely a second-order polynomial representation 

<fa(t) ~ l + ai t + a 2 t 2 (t-»0) (15) 

of the initial time series. In the language of [[32|], a\ is the "velocity" and 2a 2 is the "acceleration". The relative 
influence of the velocity and acceleration is quantified by the so-called Froude number [32] defined by 

„2 



d2 



(16) 



Working with the "direct" series (|15|), the program described in the previous section provides the sought observable, 
that we name ^(i), equal to the re-summed approximant: 



a 2 



ai 



*2(*) = $21 V) = l + ait exp 
We will also study, when necessary, the inverse function <j)2i(t) of 4>2{t), 

<j)2i{t) = 1/^2 (t) = 1 - ait + (-oa + a\) t 2 + ... = 1 + M + b 2 t 2 + 
The corresponding Froude number Fj is 



Fi 



bi 



F d 



Fi 



Ft -I 



(17) 



(18) 



(19) 



b 2 a 2 -a 2 F d - 1 ' 

In cases when better convergence can be achieved by working with the "inverse" series ©, the observable 
will be expressed through the approximant corresponding to the inverse series as follows: 

-l 



m) = [tin®] 



1 + b\t exp 



The convergence of our procedure is checked by estimating the value of the multiplier 



M 21 (i) 



1 + —t I exp 

ai 



-t 



ai 



(20) 



(21) 



The same expression gives the multiplier for the inversion function by changing all the coefficients a^'s into bj's, 
where the fcj's are related to the a^s via expression (18). 
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3.2 Negative velocity a% < (downward trend) 

3.2.1 Negative acceleration a 2 < 

Working with positive observables, this case corresponds to the possibility that the observable vanishes in finite time. 
This zero-crossing occurs at i c21 given by 

l + F d Z d exp[Z d ]=0, (22) 

where Z d = ^ t C 2\. The corresponding multiplier M2i(i C 2i) is larger than 1, signaling a possible problem with the 
convergence of the functional renormalization method. This may signal an instability in the time dynamics of the time 
series close to the zero-crossing time. 

a 2 

3.2.2 Moderate positive acceleration < a 2 < -j: 

The observable goes again to zero in finite time. The time £ c2 i at which the observable vanishes is given again by 



(22). The corresponding multiplier M\(t C 2\) is now less than 1, signaling a stable scenario. 



3.2.3 Strong positive acceleration -± < a 2 

The acceleration is sufficiently positive to counter-balance the negative trend and the observable goes through a 
minimum before rebounding upward. The time t min at which the minimum $ min = 1 — — is reached is given 

by 

^min • (23) 

(12 

This is a very stable situation since the multiplier is zero at the minimum. Our analysis suggests that the extrapolation 
to the future is the most credible for this situation in the downward trend case. 

3.3 Positive velocity ai > (upward trend) 

3.3.1 Negative acceleration: a 2 < 

The observable increases up to a maximum $ max = 1 — and then decreases after it. The time t max of the maximum 
is given by 

t --^ 

l max — 

This is a very stable situation since the multiplier M 2 i is less than 1 for arbitrary time t and is zero at the maximum. 

3.3.2 Positive acceleration 

This case requires the inversion of the time series, since the multipliers of the direct series are always larger than 1: the 
expansion is not convergent and does not supply any information on the real-axis singularities. The inversion maps 
these two cases on the cases previously considered with multipliers smaller than 1. In addition, the corresponding 
series has alternating signs, suggesting that finite-time singularity is indeed present. 

Moderate positive acceleration: < a 2 < a\ The observable increases up to a maximum and then decreases 
after it. The time t max of the maximum is given by 



bi oi 
6 2 a\ - a 2 



^ (24) 

"'max L 2 ' v 4 ^/ 
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and the value at maximum is 

$max = (l - ~fj ■ (25) 

This solution is very stable and thus very credible as the multiplier for 4>%u{t) is always smaller than 1 and vanishes 
at the maximum. The original function is obtained by the inversion = ($2i/(*))~ ( see a ^ so equation (|l8|)). 

Large positive acceleration: ^ a 2 < 02 or Fa < There is an upward finite-time singularity at t&i given by 

the zero of the inverted renormalized expansion, i.e., by the solution of 

021/ (W) =l + FiZi exp[Z 4 ] = , (26) 

2 

where Z{ = t&x 1 ■ The multiplier for 4>* 2X1 (l) ) is smaller than 1, signaling the convergence of the theoretical 
procedure for the inverse functions. Recall again that, in order to return to the original function, the inversion = 
(</>2i/(£)) 1 should be performed, leading to a divergence of the original observable at the point where its inverse 
crosses zero. 



3.4 Example 

Let us consider a function with known singular behavior and compare the results obtained by means of self-similar 



approximations with the exact values. We consider the function 1/ cos(i) which was found [33] to describe the time- 
evolution of the crack length of a self-consistent model of damage and used [|TJ] for prediction tests. Starting from 
t = 0, the function 1/ cos(i) possesses a singularity at t = tt/2 « 1.5708, and its expansion up to second order in 
powers of t 2 reads 

<fa{t) ~ 1 + 1/2 t 2 + 5/24 t 4 + ... (27) 

Note that, since 1/ cos(t) is an even function, the relevant variable for the expansion is indeed t 2 and our previous 
classification must be applied to the expansion (27) up to second order in t 2 (i.e., to fourth-order i 4 ). 



The corresponding Froude number is F^ = 1.2 and obeys the condition 1 < F^ < i.e., the previous section 
shows that one can anticipate a finite-time singularity, only on the basis of the quadratic polynomial expansion. The 
expansion of the inverse of 4>2(t) reads 



<hi{t) ~ 1-1/2^ + 1/24 t 4 + ... (28) 

As expected from our previous analysis in terms of the stability of the renormalization flow quantified by the 
multipliers, the position of the critical time is better estimated from the zero of ^2i/(*c2i) = 0> which gives i c2 i = 
1.56645, compared to the value £20 = 1-59245 suggested from the condition </>20/(*2o) = 0. This last value £20 
is nothing but the prediction of the critical time from the "bare" polynomial expansion. We refer to ( Jl~4] ) for the 
definition of 4>2ii^) an( ^ 02O/(*)- This shows that using a control parameter, which ensures that the renormalized 
expansion retrieves the polynomial expansion, improves the prediction of the critical time. 



4 Illustration of the selection of scenarios by the convergence of higher-order ap- 
proximants 

We now exploit the previous example of the function 1/ cos(t) described in section [O] to illustrate how the concept of 
convergence of the approximants developed in section II can be used for an improved determination of the singularity. 
Higher-order expansions are given by 

fal(t) ~ <hi{t) - 1/720 t 6 + (63 = -1/720) (29) 
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and 

4> 4I (t) ~ 03/ (i) + 1/40320 t 8 + (6 4 = 1/40320) . (30) 

The corresponding higher order approximants are 

m) = mi(t))~ 1 = (l + ht 2 exp( b ^t 2 exp( b ^r 2 t 2 ))) ' , (31) 



with 



6 2 



t-2 = 1 - . (32) 



and 



with 



26163 

*2(t) = (^(t))" 1 = (l + 6it 2 exp(^t 2 exp(^T 2 t 2 exp(|r3t 2 )))) ' , (33) 

r 3 = - — (246?6 4 6 2 + 56* - 126 2 6 2 - 126 1 6 3 6 2 ) . (34) 

As explained in section II, the control parameters r 2 and t 3 are determined from the condition that the expansion of 
the approximants at small t coincide with the perturbative expression of the initial polynomial representation. 

For each approximant of a given different order and with a given number of control parameters, we obtain a 
prediction for the position of the singularity obtained from the condition of zero-crossing of the inverse approximant. 
Let us compare these values between them and with the critical time derived from the corresponding initial polynomial 
fit. We find the following results: 

t c21 = 1.56645, 

t c3l = 1.55134, t c32 = 1.57067, 
too. = 1.55193, t c4 2 = 1.57048, t cA3 = 1.57079; 



corresponding to the following multipliers defined by (|12[): 

M 21 (t c21 ) = 0.6484, 

M 3 i(t c3 i) = 0.68955, M 32 (t c32 ) = 0.63708, 

M 41 (t c41 ) = 0.68743, M 42 (i c42 ) = 0.63772, M 43 (t c43 ) = 0.63663. 

The values £ c2 i, t&i, t^, £ C 4l> *c42 and t c43 should be compared with the values determined from the "bare" 
polynomials 

t 20 = 1.59 2 45, t 30 = 1.56 991, t 40 = 1.57082 . 

The convergence of this last sequence can be analyzed by looking at the corresponding sequence of multipliers 
raj(i,0) defined by ©: 

m 2 {t 2) = 0.57735, m 3 (t 03 ) = 0.63985, m 4 (i 4) = 0.63651. 

Overall these results show indeed a good convergence for i c4 i , £ c42 and £ c43 by increasing the number of control 
parameters and along the diagonal t c2 \ , i c32 , i c43 by increasing simultaneously both the order of the polynomial and 
the number of control parameters. The same observation holds for the multipliers M. In contrast, the bare polynomials 
give a slower non-monotonous convergence. This confirms the improvement brought by our scheme to obtain a better 
determination of the critical time t c . 
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5 Non-local control through the moments of the function to predict 



Section II and our subsequent tests have shown that the control parameters provide a mean to improve the extrapolation 
to the future by imposing some constraint on the approximants. Up to now, we have used the constraint that the 
approximants must retrieve the "bare" polynomial expansions at small times t. This corresponds to constraints which 
are local in time. 

It is interesting and potentially useful to investigate the possibility of using more global "non-perturbative" con- 
straints. A possible example is when, either from a priori theoretical knowledge or from experimental or empirical 
measurements, we get hold of the first j — 1 moments Hi, i = 1,2.., j — 1 of the sought function <j){t) in some interval 
[0,T]: 

m = ! • (35) 

Jo 

Our notation means that, for j = 2, we know only the zero moment ("mass" or integral <j){t) from to T), for j = 3, 
we know both the zeroth and first moment, etc... 

Endowed with this knowledge of the first j — 1 moments, we can condition the control parameters ti i T2 i ..t j -i 
demanding that the reconstructed approximants have exactly the right values of their moments: 

r «(t,Tl ,V2,...,Tj-l) f^dt = m . (36) 



/ 

Jo 
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For j = 2, we have one equation for t\ : 



/ <l>2(t,Ti)dt = HQ 

Jo 



For j = 3, we obtain two equations for t\ and t 2 : 

f 4>l(t,TiT 2 )dt = ho, [ 4>t{t,TiT 2 ) t dt = hi ■ 
Jo Jo 

For j = 4, we have three equations for n, r 2 and T3: 

rT fT rT 

/ <Pl{t, Tl T 2 , T 3 )dt = HO, / <t>*S, Tl T 2 T 3 ) tdt = Hi , / (f>*S, Tl T 2 T 3 ) t 2 dt 

Jo Jo Jo 

Based on these conditions, two different problems seem most natural. The first one is to construct an approximate 
representation of the function <fi(t) in the same interval [0, T] where moments are given or measured. In the case 
where the moments are obtained through some experimental procedure leading to some measurement errors, this first 
problem amounts to filter out the noise in the measurement interval [0, T]. The second problem that we shall address 
here consists in extrapolating to times t > T. 

Using the previous example of the function 1/ cos(t), let us consider the following case T = y/2. This value is 
"natural" as it is the root of 1 + bit 2 . Constraining the control parameters by the knowledge of the moments in the 
interval [0, y/2], we obtain the corresponding approximants. The analysis of the zero of the inverse approximants give 
the following estimations for the critical times and the corresponding multipliers 

t c21 = 1.56888, t c32 = 1.57077, t c43 = 1.570796, 

M 21 (t c21 ) = 0.64388, M 32 (t c32 ) = 0.63678, M 43 (t c43 ) = 0.63662 . 

Notice the extremely good quality of the convergence of both the critical times and of the multipliers. 

It is also possible to use an hybrid approach, where some control parameters are obtained from the agreement with 
the polynomial expansion at small time t, while the remaining ones are determined from the conditions on the known 
moments. As an illustration, we show the fifth-order approximant 

(f>l I (t,Ti,T 2 r 3 ,r 4 ) = ^1 + 6 4 t 2 exp ^nt 2 exp ^r 2 £ 2 exp ^r 3 t 2 exp(^r 4 t 2 )^^ , (37) 
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where n = 1 is conditioned by the polynomial expansion and the other control parameters should be calculated from 
the system of equations 

Jo 

/ 05/(M> r 2, r 3,T"4)^ = /Xl , 

Jo 

[ <j)* 5I (t,l,T 2 T 3 ,T 4 )t 2 dt = /X 2 . (38) 

Jo 

Note that we do not even need to know the exact values 63, 64 and 65 of the polynomial expansion since they can 
be included in the corresponding controls. This results from the fact that the constraints on the moments overwhelm 
the initial information on the coefficients of the polynomial expansion. We find t c ^ = 1.570796 in extremely good 
agreement with the exact critical time t c = tt/2 = 1.5707963. 



6 Classification and forecasting of critical times beyond mean-field 



We now use the formalism of section II and relax the condition s — > +00 on the exponent of the algebraic transfor- 
mation, which amounted to impose the mean-field value z = — 1 of the critical exponent. The control exponent s will 
be determined from the optimization of the convergence and the stability of the renormalization flow, according the 
general principles developed by Yukalov and Gluzman [23]-[|3( 



6.1 General procedure 

Consider as before an expansion of an observable <p(t) in powers of a variable t (time) given by 4>k{t) = J2t=o a n t n 
where ao = 1 without loss of generality by suitable normalization. The method of algebraic self-similar renormal- 
gives the following general recurrence formula for the approximant of order k as a function of the 
i(t) up to order k — 1: 



ization [23] 
expansion <pk 



<Pt(t) = 4>k-x{t) 



1 



a k 



-s/k 



k a k 



-s/k 



(39) 



where, in general, s = Sk(t) depends on the approximation number k and on the variable t. These approximants 
automatically agree with their corresponding polynomial expansions and the sole way to impose some control is to 
restrict s using some conditions of rather general nature such as convergence of the sequence of approximants. 

In the sequel, we assume that only the second-order expansion is available. In order to determine the critical 



exponent z, we follow Yukalov and Gluzman [25] and construct the two approximants available from the knowledge 
of the two coefficients a% and a 2 . They can be readily obtained from the general formula (39). The first order 
approximant is simply 



si 



(40) 



Representing fait) as 02 (i) 
obtain the second order approximant 



1 + a\t{l + ^ t) and applying the general formula to the expression in brackets, we 



y* 2 (t) = l + ai t 1 



«2 

a\s 2 



S3 



(41) 



Let us assume further that si = S2 = s, where s is the limiting value of the control function of the algebraic 



transformation at the critical point. It is apparent from the form of (|39j) that s plays the role of the critical index z. As 



it was explained in [25], this is justified in the vicinity of a stable fixed point. 
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The condition of maximum stability of the renormalization flow is equivalent to imposing that the difference 
_ be a minimum with respect to the set of parameters (the so-called minimal difference condition). We discuss 

below an application of the technique applied to direct second-order expansions. 

In the present work, we are interested in testing for the existence of a finite-time singularity or critical point. 

Looking for such an occurrence, we need to solve the minimal difference condition which amounts to look for the 

solutions of the two equations in terms of the two variables t c and s: 

¥*(t et s) = and %(t c , s) = . (42) 
The vanishing of given by ( ftO| ) gives 

t c = a/oi . (43) 
The second condition ^ = w i tn ( pl| ) provides an estimate s for the critical index. In terms of Froude parameter Fa 



defined by (|16[), the second condition of (|42J) can be conveniently written as 

1 + a (l-F^y S = 0. (44) 

This gives s = s(Fd) as a function of Froude number. In the cases when this equation does not have a real solution, 
we determine the control parameter s (which, we recall, is more general entity than the critical index) from the 
minimization of ^^(tc, s): 

imnjlls (l-F d i y S ^J . (45) 

The Yukalov-Gluzman technique then confronts the two approximants \t* and ^J;: a f ter tne i r difference is min- 
imized, it remains to decide which one of them is the best re-summed expression originating from the original per- 
turbative expansion. Implicit in this approach is the concept that the renormalization approach might not be fully 
convergent asymptotically but only locally. Such a decision can be made based on the analysis of the corresponding 
multipliers 

yielding 

M 1 (t,8)=(l-±) , (46) 



M 2 {t,s) = [l- — t 
ais 



1 + % l-*L t 



-i 



(47) 



The more stable solution corresponding to the smallest \Mj\ should then be selected. We find in general that ^\ 
has the smallest multiplier in the critical region, which onset is determined by the condition | Mi (t,s)\ << \ M 2 (t,s)\, 
provided that a solution to (|44|) or (45) exists. On the other hand, we find that ^2 prevails in some "pseudocritical" 



regime when the first solution becomes unstable. One can make this selection process automatic by the weighting 



procedure advocated in [ |30| ] which has also been used in [|32|]. The weighting procedure amounts to defining an 
average of the two approximants with weights inversely proportional to their multipliers. The rational for this approach 
is that the inverse of the multipliers can be shown to play a role similar to the probability that the system visits the 
dynamical state described by the corresponding approximant. The resulting function is 

, s) = n(t,s)\M 1 (t,s)\- l + n(t,s)\M 2 (t, S )\- 1 
iM 1 (t, S )|- i + iM 2 (t, S )r i 
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Since we are interested in the physically meaningful case when the critical time t c = s/a% is positive (it may be 
infinite), a modification is required when s/a\ is found negative or simply where there are no real solution to (fl2|). 



We take this situation as a signal that one should use the inverse function defined by (18) as the relevant expansion to 
obtain the most stable scenario. The sought function is then defined as the inverse of the weighted-average for inverse 
renormalized approximant. The approximants and corresponding multipliers are calculated using the parameters Oj's 
(equation (15)) changed into frj's (equation (^)). The final solution reads 



«r fc ,) . *}(*, „)-i = f^(M)|M lf (M)|-; + ^(M)|M ;/ (M)|- 1 \ " 

V IMu^^r' + lMa/Ct.a)!- 1 / 

6.2 Negative velocity ai < (downward trend) 

Using the general procedure of section ^T|, we now present the corresponding classification of the different possible 
regimes. 

6.2.1 Negative acceleration a 2 < 

In this situation, the inverse function has always a singularity. This corresponds for the direct observable function to 
vanish in finite time (critical regime I) at t c = s(Fd)/ai, where s{Fd) is the negative solution of equation (p|). Both 



approximants and \P 2 contribute to the expression ( |49| ) and ^ progressively dominates at time approaches t 



Re- 



solution ^2 contribute to the average ( |48| ) more than \t| , because |M 2 | always smaller than|Mi| . As the multi- 
plier of Vl/*(i) eventually blows up to infinity at t c , ^(i) ends by dominating the behavior of ty*(t) . 

The asymptotic behavior of the average close to t c is determined by ^ an(1 * s characterized by exponent z = 
1, notwithstanding the fact that the control exponent s is fractional. This corresponds to the situation where the 
observable goes to zero linearly in time. 

6.2.2 Positive moderate acceleration < a 2 < f^, F < F d where F = f 1 — e~ x ' e 

In this region of parameter Fq = (l — e~ 1//e ^ = 3.249 < F^, the observable still goes to zero in finite time (critical 
regime I). The time t c at which the observable vanishes is given by the same formula as in section |6.2.1[ This solution 
exists as long as there is solution s{Fd) to (|4j). When Fd becomes too small, ( |44"| ) has no more any solution and this 
corresponds to the critical regime II discussed in the next section. The boundary between these two regimes occurs at 
the Froude value Fq determined by adding the condition 

d(%(s/ ai ,s)) 



ds 



(50) 



to the general equation ^(s/ai, s) = 0. The minimum, solution of (50), is located at 

1 



In 



(51) 



and coincides with the zero of ^ 2 on ly f° r the specific value of the Froude number Fq thus determined by 

F = (l- e" 1/e ) ^ = 3.249 . (52) 



As Fd — > oo, equation (44) can be solved exactly and s(Fd — > oo) = —1, which is a mean-field value z = 1. In other 
words, given in advance a linear function <p(t) = 1 — \ai\t, our technique will indeed reconstruct it! Let us expand 
Eq.(35) around this exactly solvable limit in powers of a small parameter 1/Fd = y, 

1 + s (l - F^Y 3 ~ 1 + s + s 2 y + .... (53) 
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l + yi^4y + (y-O). (54) 



Then, 

1 

S ~ 2^ 

This expression breaks down around F d = 4, where the whole equation ( |44| ) should be considered down to Fq, where 

s(F ) = -e. 



Solution starts to contribute to the average ( |48| ) more than 'J'?;, as soon as * satisfy condition |Mi| <C I-M2I, as 



shown in Fig. 1, which represents the dependence of the two approximants \£f(i) (equation (4C)) and ^(t) (equation 
(pd])) and of their weighted average ty*(t) given by (Q. As the multiplier of eventually vanishes at t c , 
ends by dominating the behavior of and the average demonstrates critical behavior with positive fractional 

exponent z = —s(F d ). Thus, in this region of F d , as t goes to t c , we obtain a critical behavior with fractional s 
playing the role of critical index z. It means that the exponent is now different from —1 and is determined by the 
solution of eq. (|44|). 

a 2 

6.2.3 Strong positive acceleration < a 2 , F d < F 



The critical regime (I) is now replaced by the critical regime (II), such that equation ( J44| ) has no more any solution 



and the control exponent and critical time are determined from the minimization of eq. (45), which gives 



. , 1 

s = s min (F d ) = — j— — r- . (55) 

The critical index is z = —s m i n (Fd), leading to a logarithmic correction to the mean-field value. The critical time is 
given by 

, S m i n (i^) 

t c = • (56) 

ai 

Figure 2 shows the dependence of the two approximants ^*(t) (equation (^)) and ^(i) (equation (fnj)) and of their 



weighted average given by (48). The new feature is the existence of a minimum at time i m i n of *|(*) an ^> 

(57) 



therefore of a non-monotonous behavior also of the average ^f*(t), with t m j n given by 

CL\ Smin 



Trim — 

«2 Smin 



corrected by the ratio Smi " compared to the mean-field result (23) of section 3.2.3. The value of ty%(t) at t m - m is 



*5nd» = 1-^(1 + L{F d ))-^ F ^W , L(F d ) = In (-^ . 

The trajectory shown in figure 2 is rather unusual since after spending some time close to the local minimum 

of a non-critical branch ^(t)' ^ e system suddenly breaks down towards the "critical" branch ^f*(t), which then 
ends at a critical point t c . This means that the critical behavior with exponent z = —s m - m (F d ) has not disappeared 
yet. The drop occurs at a crossover time t = t cros defined as the solution to the equation \M\\ rj \M 2 \, with 
a magnitude A = ^ 2 (t C ros cros)- This regime is found for the Froude interval Fqi < F d < Fq where 

Fqi = (l — e _1 ) 1 = 1.582 is the solution of the equation 

1 + =-^=- = , (58) 

corresponding to the Froude value at which the multiplier M\{t, s) changes from stable (Mi < 1) to unstable (Mi > 
1) behavior. As the regime 1 < F d < Fqi (pseudocritical regime I) sets in, the multiplier Mi(i, s) becomes larger 
than 1, increases with time and diverges at t c . Rather than converging to the ^\(t) approximant, the weighted average 
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if?*(t) exhibits a fast change of direction to reach ^ 2 (t) at t c , as shown in figure 3. The critical branch has disappeared 
as the non-critical branch ^ 2 dominates. The presence of the approximant scenario ^| is felt only in the existence of 
some oscillations decorating the ^ 2 scenario. This regime exists for 1 < < Fqi. 



For Fd < 1 ( pseudocritical regime II), the minimum of ( ]45[ ) disappears and there is no solution either to ( |44| ) 
or to ( p3] ) anymore. This implies that we should use the inverse expansion in terms of the coefficients b\ > 0, b 2 < 
(Fi < 0). The corresponding control exponent s is then obtained from the condition 

1 + 8 (l-Ff 1 )'" = 0, (59) 

which has a negative solution s = s(Fj) leading to t c = s/b\ < 0, which is not allowed. However, similarly to 
the previous strategy to replace ([P]) by (p5|), we can look for the solution which minimizes the left-hand-side of (|59^), 
which gives s — > 00. This corresponds to a pseudocritical regime which is reminiscent of the last phase of the previous 
regime, but with formally infinite t c . In the limit s — > 00, we obtain 

*Jj(t) = exp(M), M u (t) = exp(M) 
*2j(<) = 1 + M exp(^t), M 2/ (t) = f 1 + ^ exp fe) . 

Note that ^l/W * s qualitatively similar to the mean-field solution of section III, derived for the same region of 
parameters. The multiplier Mu(t) is always larger than 1 which implies that the scenario always dominates in 
the weighted average, although the contribution of is responsible for the existence of an extra- minimum in the 



trajectory of \£* given by expression (49), as shown in figure 4. 



6.3 Positive velocity ai > (upward trend) 

6.3.1 Negative acceleration: a 2 < 

One can still define the control exponent s from the condition (0), but t c becomes negative which is undesirable. As 
in the previous section, we turn to the next possibility which is to use (|3|), whose only solution is s — > +00. This 
solution turns out to minimize the difference between two approximant scenarios. This regime is the inversion of the 
pseudocritical regime II just presented above (only with F^ < instead of Fj), as it is described by the following 
solutions, 

= exp(aii) , Mi(i) = exp(aii) . (60) 

Kfi) = 1 + ait exp(— t), M 2 (t) = ( 1 + —t) exp ( —t) . (61) 
ai V «i / \ai / 

Note that ^C^) i s nothing but the mean-field solution of section III, derived for the same parameter region. The 

multiplier M\(t) is found to be always larger than 1. Therefore, ^2 dominates in the weighted average trajectory 

*I'*(48). The contribution from ^| induces splitting of the mean-field maximum (at t = —01/02) as shown in figure 

5. 

6.3.2 Moderate positive acceleration: < a 2 < ^j+^- a\; < F d 



In this case, although equation ( |44D has a solution for s < 0, the corresponding t c is negative. After inversion of the 
initial series, this region of parameters is equivalent to 1 < Fj < F . This regime corresponds to the inverse of the 
critical regime II described above and can be described similarly. 

Consider first the region of F i < Fj < F . There is a minimum of the curve and the shape of the observable 
(\[/*) -1 ( see equation (49)) is significantly non-monotonous, due to contribution from The time t m i n of the 



minimum of (\&2*.r) 1 (maximum of ^ 2 i) i s gi yen by 

02 Smin — J- 
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where 



is the solution of the minimization 



Smin — f tp i \ (63) 



mm 



i + .(i-J7')-). 



The maximum value is (^/(min)^ w here 



*2i (mi n) = (l - Fj (1 + m )r il+m))/m) ) , = ^ (^) • (64) 

The trajectory has the same topology as shown in figure 2 except that it is the inverse of the function shown in figure 



2. The observable is predicted as the weighted average scenario given by (49) and is shown in figure 6. The "critical" 
branch shapes the weighted average ^fj as t c is approached. The weighted average scenario goes to infinity in 
finite time t c = s m i n /bi, with negative z = s m i n describing the power-low divergence. In terms of the coefficients 
dj of the polynomial expansion, this regime holds for a\ < 02 < ^^T 1 a\, i.e., for < Fj < j^tj 

(Sf^ = 0.368, = 0.692). These conditions are equivalent to F 01 < Fj < F . 

As F/ becomes smaller than Fqi, the multiplier Mu(t,s) changes from stable (My < 1) to unstable (My > 1). 
As a consequence and similarly to the change from figure 2 to figure 3, the weighted average scenario changes 
dramatically and does not exhibit anymore a critical divergence at t c . The scenario is felt only in the creation of a 
few oscillations around ^21- This regime holds for 1 < Fj < Fqi and mirrors the pseudocritical regime I previously 
described. In terms of initial coefficients, it corresponds to < a<i < F °p m 1 a\. 

6.3.3 Large positive acceleration: of < 02, Fj > Fq 

In terms of the inverse Froude number, this regime holds for Fj > Fq. The observable goes to infinity in finite time 
at a critical time t c , which is determined from the condition that the inverse quantities cross zero. The corresponding 
control exponent s(i*j) is the solution of 

1 + s (l - ~ S = , (65) 

and 

t c = s(F r )/6 1 . 

The existence of the finite-time singularity holds as long as F[ > Fq. This condition can be re-expressed in terms of 
the direct Froude number and gives F^ < -p^rr = 1-445. This regime mirrors the critical regime I. Given a second- 
order expansion <p2{t) — 1 + |di| t + |ai| 2 t 2 of a simple pole 0(f) = (1 — |ai| (with Fj — > 00), our technique 
will indeed reconstruct it! The weighted average scenario ^*j(t, s) _1 goes to infinity in finite time t c , with negative 
z = s(Fj) describing the power-low divergence. In Fig. 7 we demonstrate the dependence of the two approximants 



^l(t) and ^2(0 and °f tneir weighted average $>*(t) given by (49) . 



7 Concluding Remark 

Starting from a representation of the early time evolution of a dynamical system in terms of the polynomial expression 
of some observable (j)(t) as a function of time, we have investigated the conditions under which this early time 
dynamics may lead or may not lead to a finite-time singularity. The corresponding classification has been performed 
from the point of view of the functional renormalization method of Yukalov and Gluzman [18]-[[3(J, with the purpose 



of identifying the most stable scenarios, given the early time dynamics. The direct extension of this work is to test our 
predictions empirically, following the methodology of fli6|] developed for a particular case. 
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Figures captions 



Fig. 1: First-order approximant ty\ (equation (BO)), dashed line), second-order approximant ^ (equation (|T]), 
dotted line) and their weighted average given by ( |48| ) (continuous line) as a function of time, for positive moderate 



acceleration < a-i < Fq < where Fq = (l — e 1//e ^ 



Fig. 2: Dependence of the two approximants Vl/*(t) (equation (40), dashed line) and ^(i) (equation (|4l|), dotted 



line) and of their weighted average given by (48) (continuous line) in the regime of strong positive acceleration 



Fig. 3: Same as figure 2 in the regime 1 < Fj < Fq\ (pseudocritical regime I). Note that the weighted average 
ty*(t) exhibits a fast change of direction to reach VPjjW at t c . 



Fig. 4: Same as figure 2 for F^ < 1. 



Fig. 5: Same as figure 2 for a negative acceleration 02 < 0. 



Fig. 6: Same as figure 2 in the regime of moderate positive acceleration: < 02 < ^p-l «f; < F 



Fig. 7: Same as figure 2 in the regime ^ 1 a{ < a%, Fi > Fq. 
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